function plot_acc_sn_sp(fn_cell, legend_cell, fig_name)
n_f = length(fn_cell);
n = 100;
for i=1:n_f
    fn_cell{i}
    r{i} = importdata(fn_cell{i});
    if n > size(r{i},1), n = size(r{i},1); end
end

markers = {'-dr', '->b', '-<g', '-*k', '-sc', '-om', '-py'};

h = clf;
subplot(1, 3, 1);
for i=1:n_f
    plot(1:n, r{i}(1:n,2), markers{i}, 'LineWidth', 2); hold on   
end
ylabel('a = (sn+sp)/2');
xlabel('order')
ylim([85 95]);

subplot(1, 3, 2);
for i=1:n_f
    plot(1:n, r{i}(1:n,3), markers{i}, 'LineWidth', 2); hold on   
end
ylabel('Sensitivity');
xlabel('order')
ylim([85 95]);

subplot(1, 3, 3);
for i=1:n_f
    plot(1:n, r{i}(1:n,4), markers{i}, 'LineWidth', 2); hold on   
end
ylabel('Specificity');
xlabel('order')
ylim([85 95]);
h_legend = legend(legend_cell, 'Location', 'SouthEast')

set(gcf,'NextPlot','add');
axes;
h_title = title('Short D.melanogaster CDS (Length = 60, T = 200)');
set(gca,'Visible','off');
set(h_title,'Visible','on');
set(h_legend,'Visible','on');


print(h, '-dpdf', fig_name);
print(h, '-depsc', [fig_name '.eps']); %'drosophila-60-accuracy-prune-0.001.pdf';
end